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PROPAGATION DELAY IN THE ATMOSPHERE 

David M. LeVine 
University of Maryland 


ABSTRACT 

The structure of the lower atmosphere and the ionosphere are dis- 
cussed and an introduction to ray tracing in the atmosphere is given. 
Models for each region appropriate for satellite-to-earth communi- 
cation in the frequency range 1 GHz to 15 GHz are discussed. As- 
suming a spherically symmetric atmosphere, propagation delay is 
computed for the lower atmosphere and for the ionosphere. The 
total atmospheric delay is also computed as a function of atmospheric 
conditions, frequency and altitude of the signal source. 
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PROPAGATION DELAY IN THE ATMOSPHERE 


I. INTRODUCTION 

The problem which motivated this investigation was an analysis of the limits 
which the atmosphere places on a one-way satellite-borne system for clock 
synchronization. Presumably an atomic clock is to !>e placed aboard a satel- 
lite to beam timing information to the ground. Ideally such a clock would be able 
to beam the timing information to stations which are widely distributed over the 
globe and which arc equipped with only minimal receiving equipment. However, 
in such a one-way timing system, a correction for the finite velocity of propaga- 
tion of the signal is necessary in order to obtain practical synchronization. For 
example, in vacuum, it would take about 3 milli-seconds for a signal to reach 
the ground from a satellite 1,000 km above the receiving station; however, 
present demand exists for micro-second synchronization and present state of the 
art is about one-tenth of a micro-second. An obvious first order correction to 
account for the propagation delay is to use the straight line distance between 
satellite and receiver (slant range) divided by the speed of light in vacuum. 
Presuming that the relative position of satellite and receiver is known with 
sufficient accuracy, this provides a simple correction easily implemented by the 
receiver. A major objective of this report is to provide sufficient information 
about propagation delay in the atmosphere to permit a quantitative evaluation of 
such a correction scheme. 

This report is intended to be tutorial and to lay the basis for continued research 
as much as to be a vehicle to report on the author’s work. Consequently, the 
first three sections of this report are devoted to background material consisting 
of a brief summary of the general structure of the atmosphere and the back- 
ground for simple ray tracing. Emphasis is placed on available models of the 
atmosphere because they have been used in the past, and will continue to be used 
for predictions until an improved model, w ith easily determined parameters, is 
developed. The remainder of the report is devoted to the calculations made by 
the author. Calculations of time delay have been made using up-to-date models 
for a variety of atmospheric and geometric conditions. Data is presented for the 
propagation delay in the lower atmosphere, the ionosphere and for the composite 
atmosphere. 
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II. THE LOWER ATMOSPHERE 


A. Introduction 

The atmosphere can be separated into many regions. This is done most fre- 
quently on the basis of the temperature gradient of the atmosphere, on the basis 
of its composition, or on the basis of its electrical properties. 

Between sea level and about 100 km the temperature profile of the atmosphere 
consists of alternate regions of increasing or decreasing temperature separated 
by regions of relatively constant temperature. For example, between the earth's 
surface and about 6 to 15 km, depending on the season and location, the average 
temperature decreases at about 0.5 K per kilometer. This region is referred 
to as the troposphere and is that region of the atmosphere which is directly as- 
sociated with rain, clouds and those other phenomena generally called weather. 
Above the troposphere is a thin region of relatively constant temperature called 
the tropopause which is followed by a region of slowly increasing temperature 
called the stratosphere. The stratosphere is a region of generally laminar flow 
and is characterized bv the relatively abundant presence of ozone, and it ends 
at about 50 km in another region of relatively constant temperature called the 
stratopause. Above the stratopause is the mesosphere, a region of high velocity 
turbulent winds, in which the temperature once again decreases up to about 80 km 
where it ends in another thin region of constant temperature called the mesopause. 
Above the mesopause lies a thick region of ill defined extent called the ther- 
mosphere in which the temperature increases. However, at the heights of the 
thermosphere the composition of the atmosphere is no longer uniform, ioniza- 
tion by solar radiation becomes significant and the mean free path increases 
significantly. Consequently, care must l>e exercised when referring to "tempera- 
ture" in this region. See Figure 1 for a typical temperature profile for the at- 
mosphere^^ and definition of the regions. 

Composition lly, the atmosphere is relatively homogeneous below the mesopause 
(about 80 km) and this region is often referred to as the homosphere. The ho- 
mosphere consists, by volume, of about 7 8.08 percent of molecular nitrogen, 

20.04 percent of molecular oxy gen, .03 percent argon and traces of such elements 
as helium, krypton, neon and xenon. The mean molecular weight in the homo- 
sphere is about 28.07. Another important constituent of the atmosphere is water 
vapor; However, its content is extremely variable although the amount of water 
vapor decreases rapidly with altitude and is almost negligible above the tropo- 
pause. Above the homosphere the elements of the atmosphere begin to separate 
by diffusion and the chemical composition is changed due to dissociation caused 
by solar radiation. This region is often called the heterosphere. In the hetero- 
sphere molecular oxygen becomes dissociated by solar radiation and rapidly be- 
comes the dominant component of the atmosphere, but eventually the lighter 


elements, helium ami hydrogen (produced by dissociation of other atmospheric 
gases) diffuse upward and become dominant. In fact, evidence exists to suggest 
the existence of a region dominated by the presence of helium and lying some- 
where between 1,000 and 2,000 km and given 'he name heliosphere. The region 
above the heliosphere is predominately cotr^josed of hydrogen. At heights of 
about 700 to 1,500 km the mean free path has increased sufficiently and kinetic 
energy of the constituents has increased sufficiently that a significant portion of 
the atmospheric gas can escape the earth's gravitational field. However, the at- 
mosphere as such does not end here but rather merges gradually with the ex- 
tended solar atmosphere. The region out to many earth radii but within the 
magnetosphere contains protons and electrons emitted from the sun and this 
region is sometimes called the protonosphere. It extends down to between 1,000 
and 5,000 km. 

In terms of its effect on wave propagation the atmosphere seixirates naturally 
into two regions which are distinguished by the occurrence of photoionization of 
the atmospheric gas. Above about GO km the atmosphere is very sparse as com- 
pared to the surface and would for all practical purposes l)e negligible from a 
propagation point of view were it not for the fact that at this altitude a signif cant 
portion of the gas becomes ionized by solar radiation. The resultant plasma can 
have a great effect on wave propagation, and this region is called the ionosphere. 
The atmosphere below GO km will be called the ’’lower atmosphere” and is dis- 
tinguished from the ionosphere in that it is a relatively dense gas of neutral 
particles. The lower atmosphere contains practically all of the atmosphere, the 
density at 60 km being about 10 4 of its surface value. 


B. The Electrical Properties of the Lower Atmosphere 

The relative dielectric constant of a gas is given in terms of the polarizibility, 
« 0 , of the individual molecules by the Clausius-Mossotti relationship: 
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where is the relative dielectric constant, N () is Avogadro’s number, M is the 
molecular weight and p is the density of the gas. If, in addition, the molecules 
have a residual dipole moment, m )t the Clausius-Mossotti formula can be 
generalized to account for the alignment of those dipoles by the external field. 


4 



In this case: 
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where < is the frequency of the external field and is the relaxation time. For 
frequencies less than about 100 GHz the term -t is negligible for air.( 2 * 3 ) Ne- 
glecting the loss term and assuming that air behaves like an ideal gas (i.e., is 
proportional to the ratio of pressure to temperature) and that * r ^ 1* one finds 
that: 
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where P ( is the pc^tial pressure of each non-|)olar component of air and is the 
partial pressure of the components with polar molecules. (The only significant 
com’>onent with polar molecules is water vapor and the important non-polar com- 
ponents are "dry air" and C0 2 . The K t , A ( and are constants and T is the 
temperature. 

Equation II— 3 can be solved for the index of refraction, 77 ~t ' , and for air 

the approximation p r = 1, e f - 1 « 1 are good. Hence 77 /"7~ 1 4 (1/2) (< r - l) 

is a good approximation for air. However, rather than use the index of refraction 
directly it proves convenient to define a new quantity, the refractivity, N, by the 
relationship: N (r/ - 1) x 10 f> . Then, it follows from Equation II-3 that: 
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The constants A , B and C depend on the molecular properties of the atmospheric 
gases — on their polarizability and permanent dipole moment — and are obtained 
by measurements of the atmospheric gas. P is the total pressure of air and P w 
is the partial pressure of water vapor. The most frequently used values for the 
constants are due to Smith and Weintraub.^ Using their values: 
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whore T is measured in *K and pressure is measured in millibars. This formula 
is accurate to within l/2 percent; however, it is usually approximated with small 
a rror by: 
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Note that Equation II— 6 is frequency independent. The lower atmosphere appears 
to be independent of frequency below about 30 GHz^ when dispersive effects due 
to water vapor and oxygen absorption lines become important. 

C. Models for the louv er Atm os phere 

Ideally one can obtain the information about refractivity for propagation studies 
from Equation II-6 and measured values ofT.P and P w . Unfortunately, for 
purposes of satellite communications, knowledge of these parameters is required 
over a very great range of vertical and horizontal distances, over which simul- 
taneous measurement of these parameters is difficult. One alternative approach 
is to build an idealized model of the atmosphere starting, for example, with the 
assumption that the atmosphere is an ideal, stationary gas and using a realistic 
model ior its temperature profile and the gravitational field. With some mean 
parameters such as surface density and composition as boundary conditions, this 
model could be used to construct, analytically, profiles for other atmospheric 
parameters. This data could be used in 11-6 to obtain an "ideal” model of refrac- 
tivity as a function of altitude. One such model is the U.S. Standard Atmosphere, 
1962. Unfortunately, this is not a very satisfactory method because even the 
average diurnal and seasonal changes in atmospheric properties are not repre- 
sented. A somewhat more feasible approach is to observe N(h) at many locations 
over long periods of time and then construct models for the refractivity based on 
the long-term mean values, perhaps with a parameter or two with which to tie 
the model to geographical and temporal changes in the mean properties. A num- 
ber of such models exist. 

The earliest of these models were based on the following general characteristics 
of the refractivity: 

1 . The refractivity varies between about 200 and 450 with 313 being the 
mean value for the continental U. S. 

2. Between the surface and 1 km, the refractivity is roughly a linear, de- 
creasing function of altitude. 
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3. There is a strong correlation betwoen the mean surface value of 
refractivitv, N s , and the change in N between the surface and 1 km: 

\ < where 7.32 and . o u . , r, 7 7 lor tin !'. s. 

4. The value of N at 9 km above sea level is extraordinarily constant at 
about 105. 

Among the models for the refractivity which are based on these properties are 
the following: 

a. Linear Model: In this model N(h) N s + 2N(h- h^ ) where b is the 
height above sea level and h is the height of the surface above sea level, both 
measured in kilometers. A linear function of altitude permits use of an ’’effective 
earth's radius" model to facilitate ray tracing, and is a good approximation close 
to the earth’s surface. It is useful for line-of-sight communications between 
earth-bound stations but is not accurate enough at high altitudes to be used in 
satellite- to-earth propagation delay studies. 

b. The CR PL Reference Atmosphere (1958) Model: This model has three 
parts: 

N(h) = N s * AN(h-h s ) h # £ h < h > + 1 

. v “a T h - h “J " 1 

= (N s ♦ AN) e 1 h ( + 1 5 h < 9 

= 105 e' b ^ h ” 91 h > 9 II-7 

where a = [i/(8 -h g )J (N s + AN)/105j and b = .1424. This model has the ad- 
vantages of allowing an ’ effective earth's radius" approach for low altitude uses 
and of providing a good match to the data for higher altitudes. Its coefficients 
have been chosen from a least squares match to the Rocket Panel data.^ 6 * 

c. The CRPL Exponential Reference Model: In this case N(h) is approxi- 
mated by a single exponential: 
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This model has the advantages of relative simplicity and continuity of dN dh. 
However, it is tied to the low altitude values of rcfractivity (i.e. N a and N) and 
is not as accurate a representation of N(h)at higher altitudes as is the CRPL 
Heferencc Atmosphere.!^* 

Each of the preceding models really depend on only one parameter, the surface 
value of refract! vity, N s , because the relationship N - v v is used either 
explicitly or implicitly, This was done intentionally to tie the models as simply 
as possible to the local atmospheric conditions. For example, such models 
could be used to correct for atmospheric effects on propagation at receiving 
stations equipped with only the minimal amount of metcrological equipment 
needed to measure N s and, particularly, without detailed monitoring of the upper 
atmosphere. On the clner hand, the availability of only one parameter limits 
the accuracy with which the models can be made to fit the mean rcfractivity 
profiles, 

d. The Hi- Exponential Model: The index of refraction depends explicitly 
on the partial pressure of w r ater vapor, yet none of the preceding models ex- 
plicitly take this dependence into account. The bi-exponential model does so 
while, at the same time, increasing the number of parameters available for 
curve fitting. This model consists of one exponential for the "dry-term" in 
Equation IT-6 and one for the "wet-term." Thus: 
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where D is the rcfractivity of dry air at the earth’s surface, W is the rcfractivity 
due to water vapor at the earth's surface, H fj is the scale height of dry air and 
H w is the scale height of water vapor. 

There is some theoretical justification for this model because assuming that the 
lower atmosphere can be approximated by an ideal gas in a uniform gravitational 
field, one finds that at constant temperature the pressure varies exponentially 
with altitude. In such an atmosphere, Equation II-G has the form of the bi- 
exponential model. A difficulty with the model is that the two scale heights, H, 
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and H , are not easily obtained: D and W can !>e obtained by on-site measurements 
of pressure, temperature and humidity, but H , and H can only be obtained by 
curve fitting to detailed measurements of N(h). On the other hand, maps of 
average values of the scale heights can be prepared.^ Of course, th »3 value of 
such a system depends on their diurnal and seasonal variability. Typically 
8.0 ii j 11 and 1,6 H. I # 0 fop tin* U. * s . (S| Both 111 gone r:d I\ higher in 
summer than in winter, l) 77.G (P ,/T ) and W 3.73 x 10*(P w /T 2 ) iRIi/lOOj 
where RH is the relativity humidity in percent. (It is to be noted that some dis- 
cussion exists as to the global validity of the bi-exponential model. H)) 

e. The Compound Bi- Exponential Model: This is among the latest and per- 
haps the most representative in the sequence of models for the refractivity. The 
lower atmosphere consists of two temperature regimes: the troposphere and the 
stratosphere. The compound bi-exponential model, by adding two additional 
parameters to the bi-exponential model, attempts to take the different atmospheric 
properties of these two regions into account while at the same time keeping the 
virtues of a separate treatment for wet and dry terms. Thus, the compound bi- 
exponential model consists of: 
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where D, W and H are as defined for the bi -exponential model and H fl is the 
scale height for the troposphere, H ( ,^ is the scale height above the troposphere, 
and h t is the height above sea level of the bottom of the tropopause (top of the 
troposphere). 

The scale heights and h t must 1 x) determined from available data, for example, 
by making a least squares fit; consequently, this model is strictly applicable 
only where complete data is available. However, compilation of refractivity 
profiles over a five year period have been achieved at a number of stations and 
the 5-year means used to obtain world wide maps of the various parameters re- 
quired for use of tb compound bi-exponential model. (12) 

Because it incorporates more of the physical structure of the atmosphere and 
because it has more parameters, the compound bi-exponential model is the most; 
representative of available models. It represents, usually to within one or two 
units of refractivity (standard deviation), the observed mean profiles, N(h). The 
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representation is often better than this; however, there are a number of regions 
on the globe where the representation is much worse. These lie within a strip 
*30° latitude on either side of the equator, predominately over the seas and over 
India. < 12 ) 


in. THE IONOSPHERE 


A. Introduction 


Compared to the lower atmosphere the ionosphere is a region of relatively low 
density and heterogeneous composition. The maximum density in the ionosphere 
is about 10 g/cc and at 400 km less than 10 54 g/cc; whereas, the density of 
the earth’s atmosphere at sea level is about 10 3 g/cc. In addition, beginning 
very near the bottom of the ionosphere the atmospheric gas is no longer kept 
homogeneous by turbulent mixing but begins to separate by diffusion into layers 
of the constituents. The i>ercentage of nitrogen decreases and that of oxygen 
increases. Near the top of the ionosphere the presence of helium becomes im- 
portant, forming a layer between about 1,000 and 2,000 km; and beyond the 
heliosphere hvdrogen is the predominate element. 

The charged particles in the ionosphere are primarily the result of ionization 
by the short wavelength (less than 2,000 A) end of the solar spectrum. The major 
contributors are singly ionized atomic and molecular nitrogen and oxygen and 
nitric oxide. However, especially near the top of the ionosphere, the concentra- 
tion of ions and electrons is appreciably affected by the influx of particles from 
the sun. 

Early investigators of the ionosphere thought that it was composed of layers and 
the names D, E and F apparently applied by Appelton^ 1 **) came into wide use. 
However, what was originally mistaken for layers were in reality' relative peaks 
in the electron density profile. Even so, the ionosphere is still separated into 
a number of regions or layers associated with these peaks. They are: 

1. The D- Region: The D-region is the lowest-most region of the ionosphere, 
extending roughly from 60 to 90 km. Surprisingly enough it is a region about 
which relatively little is known because of the difficulty of making measurements. 
Typical electron densities vary from about 10 - 100 electrons/cc at 70 km to 

10 3 - 10 s electrons/cc at 90 km. Occasionally large increases in D-region ioni- 
zation occur which apparently correlate with different types of solar disturb- 
ances/ 1 "*) The D-region is the ionospheric region of highest attenuation, and the 

o 

primary sources of ionizing radiation appear to be Lyman-a (1216 A), cosmic rays 
and some x-rays (2-8 A)/ 15 ) 

2. The E- Region: The E-region lies jusc above the D-region, from about 

90 km to between 120 to 140 km, with the peak density generally occurring at about 
110 km. The E-region was the first ionospheric region to be observed (by 
Appelton and Barnett, 1925) and has been extensively studied since then. It is 
an unusually regular region for the ionosphere although there appear to be ir- 
regular increases in electron density within the E-region whose origins are not 
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well understood. These phenomena arc referred to as sporadic- K and the other- 
wise well behaved region is called the normal E-region or just the E-region. 

The normal E-region apparently responds to the soft X-ray portion of the solar 
spectrum (10 - 170 A) and to Lyman- * (102G A) which ionize Electron 

loss is by recombination, and a Chapman layer(16, 17) with recombination loss is 
a good approximation to the normal E- region. The peak electron density of the 
E-region is about 10 s electrons/cc with 50-60% variation with solar activity, 
and to a good first approximation the maximum electron density is given by:^ 8 * 14) 


N m = 10 4 ^(180 * 1.44 R) cos* 
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where R is the mean Zurich sunspot number and * is the solar zenith angle. At 
night the electron density decreases sufficiently to make the E-region difficult 
to observe. The geographical variations in the morphology of the E-region are 
slight. 

3. The F- Region: The F -region is the ionospheric region of greatest 
electron density. It extends from the top of the E-region to between 700 and 
1500 km, and it is frequently sub-divided into two regions: F, and F 2 . 

The Fj region is located at the bottom of the F-region near the E-region. Its 
peak electron density generally occurs near 160 km and it has a half-thickness 
of about 50 km. The F t region is generally well defined only in the summer 
during the day. When it is evident, it behaves very much like a Chapman layer 
'•1th electron loss via recombination; however, there are many seasonal and 
geographical anomalies. The peak electron density is given to a first approximation 
by:(18) 


N m = 1 .24 x 10 4 f4.3 + .01 R] 2 cos 4/10 x' 
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Ionized molecular oxygen appears to be the major contributor to the Fj -region. 

The F 2 -region is a region of large electron density and large variations. A peak 
density of 2.5 x 10 s electrons/cc at a height of about 250 km appears to be char- 
acteristic of the quiet F 2 -region, and under quiet conditions the F 2 -region is 
largely represented by a recombination Chapman layer. I 19 ) However, the 


F 2 -region is a region of great variability and is strongly tied to the geomagnetic 
field. In addition, it is affected by tides and is apparently subject to strong 
internal electrical forces. Strong fluctuations in F , are more often associated 
with "magnetic storms" than electromagnetic outpourings of the sun, and the peak 
electron density may be as high as 10 electrons/cc. Among the anomalous be- 
haviour of the F 2 -region is the occurrence at middle latitudes of larger electron 
density in the winter than in the summer. In addition, the phenomena of spread- F 
appears to indicate a dynamic, turbulent F 2 -region. In light of its many variations 
and anomalies, analytic prediction of the morphology of F seems impossible at 
present and mapping the only feasible approach. 

4. The Upper Ionosphere: The upper extremities of the F-region cor- 
respond to an atmospheric region of decreasing density and increasing tempera- 
ture and mean free path. If the earth were in vacuum, particles would eventually 
escape the earth’s gravitational field and the atmosphere would eventually decay 
to zero. However, the earth is imbedded in a continual stream of jxirticulate 
radiation w'hich emanates from the sun and w r hich manifests itself near the 
earth in the form of a nearly fully ionized plasma of protons and electrons, often 
called the protonosphere. 

The electron content of the ionosphere above the F , peak has been measured only 
relatively recently. Indications are that tie density at 1,000 km above the earth’s 
surface depends on the magnetic dip angle, m, and the zenith angle of the sun and 
is given to a first approximation by:^0) 
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Incoherent scatter measurements indicate that the density of electrons at 5,000 
km is of the order of 10 J electrons/cc. (20) 

Typical density profiles for the ionosphere are shown in Figure 2. 

B. Model for _the Ionosphere 

A plasma is a difficult many bodied problem because of the electrical inter- 
actions between particles. Alternate approaches to the plasma are equations 
for its distribution function or equations for the moments of the distribution 
function. Either approach is in general overwhelmingly difficult. However, in 
the special case of a zero temperature plasma and with losses proportional to 
the current density, the distribution function in velocity space becomes unusually 
amenable to analysis. In the case of an harmonic excitation whose frequency is 
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great enough to neglect the motion of the relatively massive ions in comparison 
to the electrons, one can represent the plasma by an effective conductivity 
tensor, o': 
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where Q “ co( 1 + jZ), Z ' v c /^, i^ c is the collision frequency (the loss term is ap- 
proximated by -v c J where J is the current density), a p = yq 2 n M »‘ 0 is the 
plasma frequency, u> c r |q B M is the cyclotron frequency, c t |q|(B • a )/M 
where i = x, y or z, and N is the electron density. 

The expression for the conductivity can be substituted into Maxwell’s equations 
to obtain a description for the electromagnetic fields. The resultant equations 
happen to be equivalent to those obtained for a medium with tensor dielectric 
constant T - € 0 T = e Q £l + j (ct/o* Q jj. 

The index of refraction,^, for the hypothetical dielectric medium is obtained in 
a straight forward manner by seeking solutions to Maxwell’s equations in the 
form of plane waves: E = E 0 e i(C ‘ x wt) where k • k ” k (J 2 ^ 2 . in the case that 
the wave propagates in ti:e z -direction one finds: 


= 1 - 


Y 2 

T 


(i -;z> - 2(1 -x- jz> 


y 2 

T 


III— 5 


4(1 -X- jZ) 


+ Y l 2 


where X = [co /aA 2 , Y, = (oo c /co}cos 6, Y T = (a; c /co) sin 6, and 6 is the angle be- 
tween the earth’s magnetic field and the z -axis. (2" ) This relationship is attributed 
to Appelton and Hartree and is referred to as the Appelton-Hartree dispersion 
equation. It is the result of neglecting all thermal motion of the plasma particles 
and although useful in most radio communication problems cannot be used when 
plasma modes are important. 
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for ionospheric plasmas the maximum plasma frequency is about 20 MHz and 
the maximum cyclotron frequency is about 2 MHz. The collision frequency is 
greatest in the D- region of the ionosphere, decreasing to about 5 x 10'’ at the 
bottom of the E-region. (13) In the F 2 -region the collision frequency is generally 
less than 3 x 10 * collisions per second. Consequently, at frequencies in excess 
of a few hundred MHz, the ionosphere is essentially lossless, and the index of 
refraction is given by: 
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The two solutions of this equation are called ordinary (+) and extraordinary (-) 
modes. 

In the special case of propagation along the earth's magnetic field (tt 0), and 
neglecting loss, one obtains: 
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Checking the polarization of the waves corresponding to each solution indicates 
that both the ordinary and extraordinary waves are circularly polarized but in 
opposite senses. Since the two waves propagate at different velocities, an in- 
cident plane polarized wave undergoes a change in the direction of polarization 
(Faraday rotation) in traversing the plasma. 

For frequencies above a few hundred MHz, the distinction between ordinary and 
extraordinary modes becomes quite small and is insignificant for frequencies 
above about 1,000 MHz. (22) in this case: 
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Thus, at frequencies in the GHz and higher region, the plasma is essentially an 
isotropic, dispersive dielectric. The index of refraction varies as the 
square of the electron density. 

No satisfactory model for the electron density in the ionosphere which predicts 
the diurnal, seasonal and geomagnetic variations in the density is presently 
available This is a reflection of the complexity of the problem: The ionosphere 
is a dynamic region which depends on both the electromagnetic and particulate 
solar radiation, on the tides and on internal electrical forces. However, the 
primordial forca in the ionosphere appears to be photoionization by radiation in 
the short wavelength end of the solar spectrum, and for this there is a good model 
which can be used as a starting point. 

Chapman studied photoionization in the atmosphere even before the ionosphere 
was observed. In fact, on the basis of his studies he predicted an ionosphere. 
Assuming a monochromatic source of radiation, a homogeneous atmosphere 
composed of an ideal gas with molecules of mass, M, at temperature, T, and in a 
uniform gravitational field with gravitational constant g, Chapman showed that 
the rate of photoionization in a planar earth was:(16» 1?) 
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where q () is the maximum rate of ion production, z = (h -hj/fa, h is tne altitude, 
h (i is the height of maximum rate, H kT/Mg and \ is the solar zenith angle. As- 
suming a process of electron loss such as recombination or attachment, and that 
equilibrium is reached between loss and generation, one can relate q(x» z) to 
electron density. For recombination loss one obtains: 


f 1 - z - sec ye 1 1 
N(x* z ) = N 0 exp| 2 | 


and for attachment loss: 


N(x> z) - N 0 exp { 1 - z - sec xe~ 1 } III-ll 

where N 0 is the maximum electron density. 
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That the Chapman analysis is essentially correct is reflected in the strong re- 
semblencc of many of the ionospheric regions (especially in the quiet ionosphere) 
to a Chapman layer. This is especially so in the lower ionosphere, particularly 
in the normal E-region. To a first approximation, the ionosphere is composed 
of a series of Chapman layers roughly corresponding to each of the ionospheric 
regions. On the other hand, there are very many anomalies in a Chapman model 
for the ionosphere, which reflects the complex set of ingredients which comprise 
the ionosphere. 

One approach to modeling the electron density profile is to attempt to match a 
number of Chapman functions tc observed data, perhaps with some help from the 
rule-of-thumb formulae for )>eak densities presented in the previous section. 

An a klitional exponential function is generally required for the upi>er regions 
of the ionosphere. However, the Chapman functions are rather awkward to deal 
with, and since the data anti fit are likely to be rather approximate anyway, some 
simplifying assumptions can be made. A realistic approach(22) is to use parabolas 
to model the peak density due toD,E,Fj, and the bottom of F 2 -regions (the Chapman 
function is approximately a parabola near z - 0). The top of the ionosphere be- 
tween the peak of the F 2 -region and about 1,000 km often is Chapman-like varying 
as exp [l - h - e - *'] and can lx* so modeled; an exponential can be used for the 
upper ionosphere. 
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IV. KAY TRACING 


A. Geometrical Optics 

Consider a medium whose dielectric properties are slowly varying functions of 
position. Since Maxwell's equations admit solutions in the form of plane waves 
when - and ^ are constants, it seems reasonable to attempt a solution in the case 
of slowly varying dielectric properties with functions which resemble plane waves. 
Thus guess solutions in the form: 
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where k = a/e , c = l/^TT . Substituting Equations IV-1 into Maxwell's equa- 
tions for sourceless dielectrics, one obtains: 


Vs X F - y. c k 



X 


F 


IV-2a 


VS X H + f c F 


) 



a 


IV- 2b 


F • VS 


j “ • 7-Cnc +V • '] 


IV-2c 


H • VS 



V 'Cn y + V • H ] 


IV- 2d 


Introducing the variables J. = kx ( and the new quantities e = e Q h r y 0 
S r ks, one obtains: 
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c* • VS 
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where the operator V kV and r and h have the same dimensions, /x f and < r are 
dimensionless. 

Consequently, if the change in »* and h are small in one (free space) wavelength, 
then approximately: 
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These equations describe a limiting case referred to as "geometrical optics". 
Solving Equations IV-4 for VS, yields: 


(VS ) 2 = (VS ) 2 = t ) 2 
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This equation is referred to as the equation for the eikonal and its solution, S(r ), 
represents a collection of surfaces of constant phase. The locus of a unit vector 
perpendicular to these surfaces is called a ray. It is easily shown that for an 
isotropic, lossless medium, the time average Poynting vector is perpendicular 
to the surfaces of constant phase. Consequently, to the extent that the local 
Poynting vector represents the direction of energy flow, the tangent to a ray 
points in the direction of energy flow. 

A useful differential equation for the rays follows from their definition. (23) 

Let s be the length of arc measured along the ray path and r(s) be the position 
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vector to a point on the ray. Then dr ds is a unit vector tangent to the ray at i, 
and from the definition of a ray: 
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Therefore: 
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Differentiating both sides with respect to s one obtains the following vector dif 
ferential equation for the ray, r(s): 
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The velocity with which a signal travels along a ray — i.e. the velocity with 
which the surfaces of constant phase expand in space — can be found by differ- 
entiating the equation for constant phase with respect to time. That i3: 
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Using (d/dt ) S(r) Vs* (dr/dt) and VS(r)| " 17 one obtains: 


VS dr c 

y§| HT 7 ) V P IV-10 

where v p is the component of "phase velocity" along a ray. 

If the medium is dispersive, the velocity of propagation of a signal loses pre- 
cise meaning because the signal changes as it propagates. However, if the signal 
does not change too much, it can often be represented by a (spreading) envelope 
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which propagates at velocity: 
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where is the "center frequency" of the spectrum of the signal, is called 
the group velocity .(25) In the special case of ono dimension, one obtains: 
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Equation IV-12 should also be a reasonable approximation for slowly varying 
media in which the solutions are locally well approximated by plane waves. 

B. The i Applicability of Geom etrical Optics to the Atmosphere 

The crucial assumptions underlying the development of geometrical optics are: 

1) That the dielectric properties of the modium vary only slightly in one wave- 
length and 2) That the amplitude of the signal does not change significantly 
in one wavelength. 

The latter assumption requires the exclusion of regions where many rays cross, 
such as at focal points, and the exclusion of regions in which diffraction is im- 
portant, such as near the edges of a shadow. Neither of these regions occur in 
the spherically symmetric, slowly varying atmospheres considered in this re- 
port. The model atmospheres considered have no discontinuities in index of re- 
fraction. Also, in order that amplitude changes be small, regions of large at- 
tenuation must be neglected. However, at the frequencies considered here (1 to 
15 GHz), attenuation per wavelength in the atmosphere is small ;(24) i n fact, all 
models of the atmosphere considered in this report are lossless. 

The scale of variation of the atmospheric properties is a more complex question 
because both regions of the atmosphere have, in addition to gross features whose 
scale of variations is kilometers, local variations whose scale can be quite small. 
For example, the average index of refraction in the lower atmosphere is reason- 
ably well modeled by the compound bi -exponential model whose scale heights 
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are of the order of kilometers. Hut this model describes only the gross features 
of the lower atmosphere and is based on information which smooths out any local 
fluctuations due to rain, clouds or turbulence. Similarly, gross features of the 
ionosphere are to a crude first approximation represented by a number of parabo- 
las whose half-widths* are of the order of many kilometers. Hut on a small 
scale the ionosphere is a turbulent, inhomogeneous region, as such phenomena 
as spread- F make evident. 

The local fluctuations in cither atmospheric region are difficult to model and 
are complex phenomena whose effect on propagation is currently an active area 
of investigation. Consequently, it seems reasonable, esi>ecially in a first approach 
to the problem, to neglect the local fluctuations and treat only the gross features 
of the atmosphere. This is done here in the hopes that a discussion of the effects 
of local fluctuations will be the subject of further work. 

C. Equations for the Time Delay 

The time required for a signal to propagate between two points, Pj and p , is 
given by: 


At 



rv-i3 


where v is the appropriate signal velocity and the integration is taken along the 
ray path. The lower atmosphere is non-dispersive in the frequency regime of 
interest, so that v is the phase velocity. On the other hand, the ionosphere is 
dispersive and for frequencies above about one GHz, where the effect of the 
earth ’s magneti c field may be neglected, one has to a good first approximation, 
t] = ^ 1 - [ou p /oip". In this case, the group velocity is r 017. 

In order to determine the ray path, the assumption is made that the atmosphere 
is spherically symmetric. This first approximation means that horizontal in- 
homogeneities in the lower atmosphere such as accompany weather fronts and 
horizontal inhomogeneities in the ionosphere such as are the result of the geo- 
magnetic dependence of the F 2 -region and changes in the local solar zenith angle, 
will be neglected. This assumption has been made to permit a tractable solution 
to the problem of time delay. It should be a reasonable first approximation, 
especially if the maximum lateral dimension of the atmosphere through which 
propagation takes place is small compared to the scale of horizontal 
inhomogeneities . 


*Half the distance between the intercepts on the axis corresponding to zero electron density. 
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In the case of a spherically symmetric atmosphere, it is easily shown that the 
vector r x T 7 (r)s is constant along; the ray, (s is a unit vector tangent to the ray 
at s.) Consequently, the rays lie in planes through the origin and 7?(r)r sin? 
is equal to a constant along the ray, where ; is the angle between r and s. This 
relationship is sometimes called Bourguer's rule. 

From Figure 3 it is apparent that sin<£ - rd#/di t [r 2 ♦ (dr/d' ) 2 ]' 1 2 . Con- 
sequently, A. r rj(r) r sin i - r[r 2 ♦ (dr/d#) 2 ]" 12 where A is the constant in 
Bourguer's formula. Solving for dr/d#, one obtains: 
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Thus; 
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The unknown constant A can be expressed in terms of the initial poini ( r x ) and 
the fina’ ooint (r 2 , # 2 ) of the ray by: 
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Expressions for th< time required for a ray to travel a given path can now be 
obtained in the for t \ of two coupled integral equations. In the case of the lower 
atmosphere, using i y e phase velocity, one obtains: 
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where A is found from Equation rv- 16 . In the case of ionosphere, neglecting the 
effect of the geomagnetic field and losses, one obtains: 


At 




r dr 

|/r 2 T7 2 ( r ) - A 2 


IV-18 


where t ) /l " (^/w ) 2 and A is found from Equation IV- 16. 

The preceding two pair of equations have been solved for a number of combina- 
tions of 77 ( r) and geometry and the results are reported in the following section. 
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V. 


RESULTS 


A. The Ixnver Atmosphere 

The compound bi-exponential model was used to compute the time required for 
propagation in the lower atmosphere. The signal was assumed to start at various 
distances above the earth's surface ranging from an altitude of 50 km to 100 km, 
and for each altitude the time required for the signal to propagate to the earth 
was computed for many different paths. For each path, the actual aspect angle 
and the difference between the actual time required for propagation and the time 
which would be required in the absence of an atmosphere (free space jxith) were 
computed. The actual aspect angle is the angle which the straight line between 
sender and receiver makes with the vertical, measured at the receiving station. 
The difference between actual propagation delay and the free space propagation 
delay, which will be referred to henceforth as the "corrected time delay", is a 
significant characteristic of the path in two important respects: 1) It is an indi- 
cation of the effect of the medium on propagation; 2) It is an indication of the 
efficacy of a simple first order correction to account for propagation delay in 
timing systems. (In order to obtain micro-second time accuracy in a one-way 
timing system, some correction for the propagation delay is necessary. A 
simple first order correction is to assume that the signal propagates along the 
straight line between sender and receiver at the speed of light in vacuum — i.e. 
neglect the medium. The corrected time delay is the error incurred in such a 
correction scheme.) 

The compound bi-exponential model was used to describe the lower atmosphere 
because it is the latest and most detailed model of the lower atmosphere generally 
available at present. It is based on the five-year mean value of refractivity ob- 
served at numerous locations distributed over the globe. At each station the 
mean refractivity as a function of height has been fitted to two exponentials for 
each of four months: February, May, August and November. This data was used 
to compute time delay as a function of season for signals originating at various 
altitudes at numerous locations. Representative data is shown in Figures 4-6. 


Two curves are plotted in each figure, one for the maximum corrected time 
delay at the given location and one for the minimum corrected time delay. In 
Figure 4, the corrected time delay is plotted as a function of actual aspect for 
the mean atmosphere above Barrow, Alaska. Figure 5 represents the same data 
for the mean atmosphere above Clark Field, Phillipines. Figure 6 represents 
the corrected time delay versus actual aspect angle based on the mean atmos- 
phere above Puerto Montt, Chile. Each curve has been drawn for a signal which 
starts at an altitude of 100 km. There is very little change with altitude between 
50 and 100 km, and since virtually all of the lower atmosphere lies below 100 km, 
those signals which originate at 100 km are indicative of the maximum effect of 
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the atmosphere. Note in particular, that all of these curves have essentially the 
same shape, that the spread between maximum delay and minimum delay is 
relatively small, and that there is very little quantitative difference in the delay 
as a function of actual aspect angle measured at any of these locations. 

The locations used in Figures 4-6 were chosen because they are representative 
of the different climatic conditions around the world. 

By way of comparison, the corrected time delay was also computed using the 
CRPL exponential reference atmosphere. This data is presented in Figure 7, 
again for an altitude of 100 km. The three curves in Figure 7 represent the 
minimum delay, maximum delay and an intermediate case as predicted by this 
model. Note, that the shape of these curves is quite similar to that for the 
compound bi-exponential model, although, the corrected delay at a given aspect 
angle is, for each case, lower than the delay predicted by the compound bi- 
exponential model. For example, at vertical incidence the compound bi- 
exponential model predicts corrected delays between about 7.5 and 8.7 nano- 
seconds whereas the extremes for the CRPL model are about 6.2 and 7.3 nano- 
seconds. However, in the case N = 313, C e - 0.1439 (km) 1 the CRPL and 
compound bi-exponential models are quite close. 


The Ionosphere 


For purposes of these calculations the ionosphere has been approximated by a 
three part model. Below the peak electron density, it is assumed that the elec- 
tron density profile is parabolic. Between the height of the peak electron density 
and 1,000 km, a Chapman-like function, exp -4- [l - z - ze“ h ], is assumed for the 
electron density profile. It is assumed that the density at 1,000 km is given by 
N(1,000) = 8 x 10 3 [l + cos 2 m cos 2 (*/2)] where m is the magnetic dip angle and 
X is the solar zenith angle, and above 1,000 km the electron density is assumed 
to decay exponentially to the constant value of 10 3 electrons/cc at 5,000 km. 

The bottom of the ionosphere is assumed to occur at 100 km and calculations of 
the corrected delay have been made for signals which originate between 1,000 km 
and 5,000 km and are received at 100 km. 


In Figures 8-10 are shown curves of corrected delay versus frequency for a 
number of different values of peak electron density, representative of the actual 
ionosphere. Figure 8 displays data at vertical incidence, Figure 9 at an actual 
aspect angle of 30 degrees and Figure 10 at an actual aspect angle of 70 degrees. 
In Figure 11, the corrected time delay versus actual aspect angle is plotted for 
various values of peak electron density and a frequency of 1 GHz. A similar 
curve at a frequency of 8 GHz is shown in Figure 12. In each case m = x ~ 0, 

H m = 350 km and the half-thickness of the parabolic region is 200 km. 
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The variation of the data with m and x is not too groat, especially for large values 
of peak electron density. In order to illustrate this, the data in Figures 13 and 
14 was prepared in the extreme case, m 90 . In Figure 13, corrected time de- 
lay is plotted as a function of frequency for normal incidence (an actual aspect 
angle of zero degrees), and in Figure 14, the corrected time delay is plotted as 
a function of actual aspect angle for a frequency of 4 GHz. A corresponding 
curve for m = 0 degrees is also indicated on Figure 14 at the peak density of 
6.4 x 10 6 electrons/cc. In each case m 90°, H 350 km, the ha If- thickness 
of the parabolic region is 200 km and the signal originates at an altitude of 
5,000 km. (The change with altitude is similar to that indicated in Figures 8-10.) 

Other models for the ionosphere have also been used, although it is felt that the one 
used here is more representative of the actual ionosphere than any others tried. 
For example, a parabolic ionosphere yielded the data shown for vertical indi- 
dence in Figure 15. 

C. The Composit e Atmo sphere 

The data for the lower atmosphere and the ionosphere has been put together to 
predict the total propagation delay in the atmosphere. Doing so is somewhat dif- 
ficult because of the many parameters involved. For example, one must specify 
frequency and electron density profile for the ionosphere and location and sea- 
son for the lower atmosphere. However, for the purposes of predicting the total 
corrected time delay, the atmosphere changes very little from site to site and 
from month to month (see Figures 4-6). The atmosphere above Barrow, Alaska 
in May is reasonably representative of the typical atmospheric contribution to 
delay. Consequently, the calculations were made using a lower atmosphere 
modeled after that over Barrow, Alaska in May in conjunction with the model 
ionosphere described in the previous section.* The results are shown in Fig- 
ures 16-18, and in each figure the corrected time delay is plotted against solid 
angle. (Since the propagation problem is azimuthally symmetric, the locus of 
all stations which receive with the same corrected delay is a circle on the earth's 
surface. This circle is centered about the line which joins the center of the 
earth and the satellite and subtends a solid angle, 0, on the earth's surface which 
is plotted against the corrected time delay to obtain Figures 16-18.) For a given 
time delay, At, the corresponding solid angle obtained from Figures 16-18 identi- 
fies that portion of the earth's surface with which communications can be made 
at a corrected delay of less than or equal, to At . 


*ln each of the cases reported here m - X ~ 0, H = 350 km and the parabolic half-thickness is 
200 km. 
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In Figure 16 and 17, corrected delay is plotted against solid angle for a satellite 
altitude of 5,000 km for various values of peak electron density. The propagation 
frequency in Figure 16 is 2 GHz and in Figure 17 it is 8 GHz. 

In Figure 18 corrected delay is plotted against 3olid angle for a satellite altitude 
of 1,000 km, a frequency of 8 GHz and for various values of peak electron density. 

Comparison of Figures 17 and 18 indicates, as one would suspect, that for a 
given permissible time delay and frequency, the higher of the two satellites can 
communicate with more of the earth's surface. This is true because the signal 
is more nearly vertical for the higher satellite and because the electron density 
is relatively small above 1,000 km. 

The solid angle subtended by a circle of 3,000 miles diameter is indicated on 
each of the figures. Such a circle corresponds roughly to that portion of the 
earth's surface occupied by the continental U. S. 


♦ 
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VI. CONCLUSIONS 




The accuracy of the numerical data presented in this report is almost entirely 
determined by the models used for the atmosphere. In the case of the lower 
atmosphere the situation is reasonably well in hand. The lower atmosphere has 
been studied extensively, over long periods of time and most of its irregularities 
are restricted to a relatively small, near the surface, region. On an average 
basis, the compound bi -exponential model represents the local rofractivity pro- 
file to within a few percent over most of the globe. However, the model for the 
ionosphere is much less reliable. The dynamics of the ionosphere are quite 
complex and only relatively recently have the upper extremities of the ionosphere 
been available for study. Furthermore, the ionosphere encompasses a much 
more vast region of space than does the lower atmosphere. A realistic, quantita- 
tive, evaluation of the model used herein for the ionosphere is difficult; however, 
it was based on relatively recent observations of the density profiles. 

Of course, the models are only intended as first order representations of the 
atmosphere and are not intended to describe inhomogoneities and time variations 
in the media. For the most part these phenomena are second order effects; 
however, a few such phenomena should command attention. For example, the 
effect of rain and fog on propagation ielay as opposed to attentuation, and the 
e f fect of time varying, small scale inhomogeneities on the timing signal. When 
piopagation occurs over large distances, one should also be concerned with the 
effect of large scale inhomogeneities on the timing signal, as for example, 
weather fronts, and the changes in electron density which accompany the chang- 
ing, local, solar zenith angle. 

It also should be noted that in this investigation attention has been focused on a 
frequency range from about 1 to 15 GHz. This is the ’’window” between the 
absorption peak of water vapor at about 22.5 GHz and the realm in which the 
ionosphere must be treated as an anisotropic (eventually opaque) medium. At 
the high frequency end of this spectrum some attention should be devoted to 
attenuation due to scattering from rain, fog, etc., and at the low end of this 
spectrum inhomogeneities in the ionosphere may require attention. The possible 
use of higher frequencies, perhaps optical frequencies, in a one-way timing 
system may warrant some future extension of this frequency spectrum. 

In total, this report has aimed at summarizing the nature of the atmosphere, at 
outlining the techniques used to calculate propagation time delay through the 
atmosphere, and at making a reliable, but first order, calculation of the effect 
of the atmosphere on the propagation of timing signals at frequencies between 
about 1 GHz and 15 GHz. 
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Figure 9. Corrected Time Delay vs Frequency 
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